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Abstract: This paper addresses the stabilization issue for fractional order switching systems. Common 
Lyapunov method is generalized for fractional order systems and frequency domain stability equivalent 
to this method is proposed to prove the quadratic stability. Some examples are given to show the 
applicability and effectiveness of the proposed theory. 
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1. INTRODUCTION 

The past decade has witnessed an enormous interest in switched 
systems whose behaviour can be described mathematically us- 
ing a mixture of logic based switching and difference/differential 
equations. By a switched system we mean a hybrid dynamical 
system consisting of a family of continuous-time subsystems 
and a rule that orchestrates the switching among them (Liber- 
zon [2003], Daafouz et al. [2002]). A primary motivation for 
studying such systems came partly from the fact that switched 
systems and switched multi-controller systems have numerous 
applications in control of mechanical systems, process control, 
automotive industry, power systems, traffic control, and so on. 
In addition, there exists a large class of nonlinear systems which 
can be stabilized by switching control schemes, but cannot be 
stabilized by any continuous static state feedback control law 
Lin and Antsaklis [2009]. 

Recent efforts in switched system research typically focus on 
the analysis of dynamic behaviors, such as stability, control- 
lability and observability, and aim to design controllers with 
guaranteed stability and optimized performance (refer to Lin 
and Antsaklis [2009], Shorten et al. [2007] for a survey in recent 
results in the field). To be more precise, the study of the stability 
issues of switched systems gives rise to a number of interesting 
and challenging mathematical problems, which have been of 
increasing interest in the recent decade. 

Typically, the approach adopted to analyze these systems is 
to employ theories that have been developed for differential 
equations. To this respect, most results are based on Lyapunov's 
stability theory which has played a dominant role in the analysis 
of dynamical systems for more than a century. Existence of 
quadratic Lyapunov functions for each of the constituent LTI 
systems is not sufficient for the stability of switched systems. 
However, it is well known that the switched system is stable 
if there exists some common Lyapunov function that satisfies 
the conditions of the Lyapunov theory simultaneously for all 
constituent subsystems (see e.g. Liberzon [2003], Narendra and 
Balakrishnan [1994], Mori et al. [1998], Shim et al. [1998]. Al- 
though Molchanov and Pyatnitskii [1989] established a number 
of converse theorems, showing that such common Lyapunov 
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function always exists when the switched linear system is sta- 
ble for arbitrary switching, general conditions for determining 
the existence of a common Lyapunov function for switched 
systems are unknown. Likewise, a frequency domain method 
equivalent to the common Lyapunov one may make the control 
and stability analysis easier. For example, Kunze and Karimi 
[201 1] propose a frequency domain equivalent of common Lya- 
punov function based on strictly positive realness (SPR) of the 
system in order to analyze the quadratic stability of switching 
systems. 

Given this context, the contribution of our work is to bring 
together theories from several areas of control and to present 
stability issues in a unified manner for fractional order switch- 
ing systems. 

The remainder of this paper is organized as follows. Section 2 
provides a collection of important issues concerning stability 
of switched systems. The main contribution of this paper is 
presented in Section 3, i.e., the stability theory developed for 
fractional order switched systems. Section 4 gives some exam- 
ples to show the applicability and goodness of the proposed 
stability issues. Section 5 draws the concluding remarks. 

2. PRELIMINARIES 

When a system becomes unstable, the output of the system goes 
to infinity (or negative infinity), which often poses a security 
problem in the immediate vicinity. Also, systems which become 
unstable often incur a certain amount of physical damage, 
which can become costly. For the sake of clarity, a collection 
of important issues concerning stability of switched systems is 
given in this section, mainly using Lyapunov theory. 

2.1 Stability theorems and basic definitions 

The idea behind Lyapunov's stability theory is as follows: 
assume there exists a positive definite function with a unique 
minimum at the equilibrium. One can think of such a function 
as a generalized description of the energy of the system. If we 
perturb the state from its equilibrium, the energy will initially 
rise. If the energy of the system constantly decreases along 
the solution of the autonomous system, it will eventually bring 
the state back to the equilibrium. Such functions are called 
Lyapunov functions. While Lyapunov theorems generalize to 



nonlinear systems and locally stable equilibria we shall only 
state them in the form applicable to our system class. Consider 
an autonomous nonlinear dynamical system 

x(t)=f(x{t)),x{0)=x , (1) 
where x(t) e^Cl" denotes the system state vector, 3 an 
open set containing the origin, and / : 3> — > M" continuous on 
2>. Suppose / has an equilibrium; without loss of generality, 
we may assume that it is at origin. Then, Lyapunov stability 
for continuos systems can be summarized in the following 
theorems. 

Theorem 1. Let x = be an equilibrium point of (1). Assume 
that there exists an open set 3> with € 3 and a continuously 
differentiable function V : 3 — > M. such that: 

(1) V(0)=0, 

(2) V (x) > for all x G 3\{0}, and 

(3) H (*)/(*) < for all iei 

then x — is a stable equilibrium point of (1). 

Theorem 2. If, in addition, f£ (*)/(*) < for all x € ^\{0}, 
then x = is an asymptotically stable equilibrium point. 

Definition 1. (Quadratic Stability). A linear system 

x = Ax, (2) 
is said to be quadratically stable in R if there exists a positive 
definite matrix P eW x " such that, 

A T P + PA < 0. 

Definition 2. (t~ a Stability). The trajectory x(t) = of the sys- 
tem d jm 1 = f{t,x{t)) is t~ a asymptotically stable if the uni- 
form asymptotic stability condition is met and if there is a 
positive real a such that : 

V||jc(f)|| ,t < to 3 N(x(t),t <t ), h (x(t) ,t < t Q ) such that 
Vf <t Q , \\x(t)\\ < N(t-t x y a . 

t~ a stability will thus be used to refer to the asymptotic stability 
of fractional systems. The fact that the components of the state 
x(t) decay slowly towards following t~ a leads to fractional 
systems sometimes being treated as long memory systems. 

Let us consider a fractional order linear time invariant (FO-LTI) 
system as: 

D a x=Ax,xeR n (3) 
where a is the fractional order. 

Theorem 3. (Moze et al. [2007]). A fractional system given by 
(3) with order a, 1 < a < 2, is t ~ a asymptotically stable if and 
only if there exists a matrix P = P T > 0, P e K nx ", such that 

" (A r P + PA) sin (0) (A T P- PA) cos (0)1 ^ 
_(-A T P + PA) cos (0) (A T P + PA) sin (0)J ' k ' 

where <j) = ™ . 

Theorem 4. (Moze et al. [2007]). A fractional order system given 
by (3) with order a, < a < 1, is t~ a asymptotically stable if 
and only if there exists a positive definite matrix Pel" such 
that 

(-(-A)^s) r p + p(-(-A)^) <0. (5) 

2.2 Common quadratic Lyapunov functions 

Consider a switched system as follows: 

x = Ax,A e co{A\, ...,A£\ , (6) 



where "co" denotes the convex combination and A,,;' = l,...,L 
is the switching subsystem. According to Pardalos and Rosen 
[1987], (6) can be alternatively written as: 

L L 

x=Ax,A = Y J hA u VA,- > 0, £ Xi = 1 . (7) 

i=\ i=\ 

Theorem 5. (Boyd et al. [1994]). A system given by (7) is 
quadratically stable if and only if there exists a matrix P = P T > 
0, P e W x n, such that 

AjP + PAi<0,Vi = l,...,L. 



2.3 Quadratic stability in frequency domain 

Kunze and Karimi [2011] propose an equivalent to common 
Lyapunov stability conditions in frequency domain. The rela- 
tion between SPRness and the quadratic stability can be stated 
in the following theorem. For further information about the 
specification of state space system, refer to Section 3. 

Theorem 6. (Kunze and Karimi [2011]). Considers (s) andc2(s), 
two stable polynomials of order n, corresponding to the systems 
x=A\x andi = A2X, respectively, then the following statements 
are equivalent: 

(1) g)and^gareSPR. 

(2) |arg(ci(7fi)))-arg(c 2 0'©))| < f V co. 

(3) Ai and A2 are quadratically stable, which means that BP = 
P T >0eR nx " such that A[P + PA! < , AfP+PAj < 0. 



3. QUADRATIC STABILITY OF FRACTIONAL ORDER 
SWITCHING SYSTEMS 

This section will study two ways to obtain the quadratic stabil- 
ity of fractional order switching systems generalizing common 
Lyapunov functions for fractional order switching systems and 
obtaining an equivalent in frequency domain, respectively. 



3.1 Common quadratic Lyapunov functions of fractional order 
system 

Let us consider a fractional order switched system as: 

D a x=Ax,Aeco{A u ...,A L }, (8) 
where a is the fractional order. 

Theorem 7. A fractional system described by (8) with order a, 
1 < a < 2, is quadratically stable if and only if there exists a 
matrix P = P T > 0, P e W xn , such that 

" (AfP + PAi) sin (0) (AjP — PAi) cos (0)1 
(-AfP + PAi) cos (0) (AfP + PAt) sin (0) J ' 

Vi = l,...,Z,. (9) 

Proof!. System (8) can be rewritten as: 

L L 

D a x=Ax,A = £A,-A,-,VAi > 0,£A,- = 1. (10) 
Then, from Theorem 3, and (8), we have 



r P + P^#)sin(0) {J( 7 P - PJt) cos(<j>) 
(-^ r P + P^)cos(0) \jf T P + PJf)sm(<l>)_ 

L 

VA; >0,^A;= 1 



1=1 



1=1 



(AfP + PAi) sin(0) (Af P - PA,) cos(» 
(-AfP + PAi) cos(^) (AfP + PA,) sin(0) 



va>o,EA; = 1. 

i=l 

where ^# = J^lj A;A ( - and <j) = Therefore, it is obvious that 
(8) is quadratically stable if and only if 



<0, 



(AfP + PAi) sin(0) (AfP - PA t ) cos(0)" 
(-AfP + PAi) cos(0) (A^P + PAi) sin(</>)_ 

Vi = l,...,L. (11) 



Theorem 8. A fractional system given by (8) with order a, 
< a < 1, is quadratically stable if and only if there exists 
a matrix P = P T > 0, P e M" xn , such that 

(-(-A ! )^) r p + p(-(-A i )^) <0,V/=1,...,L. (12) 



Proof 2. Assume 



1=0 



= 0, the fractional order sys- 



tem (8) with order a, < a < 1, can be replaced by the 
following integer order system Moze et al. [2007]: 

i=A fZ ,A f GCo{A fl ,...,A fL } (13) 

z = C/x, (14) 



where A 



o 

A/a 



A, 1/a ' 




A) la 



and C 



/ 



[0 



11 



Writing (13) in an alternative way yields: 

L L 

z = A fZ ,Af = E VI/;. , VA,- > 0, E A- = 1 . (15) 

i=l i=l 

Therefore, assuming a positive definite matrix 2? > with 
proper size and, based on LMI method, the system (8) with 
order a, < a < 1, is quadratically stable if: 

A T f &> + &>A f <0=> (16) 

L 

Y,h(A T f .^ + &A fi )<0^ (17) 
i=i 

Aj.& , + &>A f .<0,Vi=l,...,L. (18) 

Then, it is obvious that expression (18) is satisfied if and only 
if (Moze et al. [2007]) 



(A\ la ) T P + PA\ la <0,V/= l,...,L, 



(19) 



where P is a positive definite matrix. In Moze et al. [2007] 
it is shown that condition (19) is sufficient but not necessary 
to guarantee quadratic stability. The necessary and sufficient 
condition for fractional order system is given by Theorem 4. 
Therefore, the necessary and sufficient condition for fractional 
order system is 

(-(-A,) z ^) r p + p(-(-A I ) 2 ^) <0,Vr=l,...,L. (20) 



3.2 Frequency domain stability 

In this section, a link between quadratic stability using Lya- 
punov theory and SPR properties will be provided, i.e., a con- 
nection between time domain and frequency domain conditions 
in order to obtain quadratic stability of fractional order switch- 
ing systems. 

Consider a stable pseudo-polynomial of order na as: 

d{s) = s na + 4V"- 1)a + • • ■ +dis a + d , (21) 
which corresponds to the fractional order system D a x = Ax. 
Furthermore, consider a polynomial of order In as: 



and 



A = 





■+cis + CQ, (22) 


3 tO X 


= Ax. Assign 








c = 


[C2«-l, C\ 


co] 




(23) 


D = 


[d n -i, d\ 


do] 




(24) 


~-c n - 
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-co" 




1 


••• 
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(25) 
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1 ••• 


















(26) 





••• 
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In the following, the necessary and sufficient condition for the 
quadratic stability of fractional order switching systems will be 
given. 

Theorem 9. Consider d\{s) and d2(s), two stable pseudo- 
polynomials of order n, corresponding to the systems D a x = 
A{X and D a x = A2X with order a, 1 < a < 2, respectively, then 
the following statements are equivalent: 



(1) 



an. 



arg ( det((Af - co z I) -2{jo))A x sin -^-)) 



arg (det((A^ - co 2 I) - 2(jco)A 2 sin -|-)) 



n w 
<2,V». 



(2) Ai and A2 are quadratically stable meaning that: BP = 
P T > € W xn such that 



(AfP + PAi) sin(0) (AfP- PA t ) cos(0)' 



<0, 



(-AfP + PAi) cos(0) (AfP + PAi) sm(<j>)_ 

Vi= 1,2. 

Proof 3. Consider c\ (s) and c 2 (s) are characteristic polynomi- 
als corresponding to x = A\x and x = A 2 x, respectively, where 

= 1,2. According to Theorem 



A,sin(0) A ; cos(</>) 
-A|-cos(0) A,-sin(0) 



6, the following statements are equivalents: 



(1) 



c 2 (s) 



and ^ are SPR, where c,-(s) = det(s/- A,-), i = 1, 2. 

(2) |£u-g(ci(7©))-arg(c2(j©))| < f V ft). 

(3) Ai and A 2 are quadratically stable meaning that: 3& = 
&> T > e R 2nx2n such that A\&> + ^Ai < , Af ^ + 
3?>A2 < 0. 



Now, consider di(s) and c/ 2 (s) are characteristic pseudo- 
polynomials corresponding to the fractional order systems 
D a x — A\x and D a x = A 2 x with order a, 1 < a < 2, respec- 
tively. The relation between A,- and A,- is given by (25). From 
(22), we have 

|arg(ci(;0)))-arg(c 2 (;co))| = 



n 



I arg {jcol - A , ) - arg {jtol - A 2 
& I arg (det((A 2 - a) 2 /) - 2(j'ft))Ai sin(0)) - 
arg(det((A^-to 2 /)-2(;a))A 2 sin((/)))| < |,Vft>. 
where / is the identity matrix with the proper size. 



Define, & 



A 



P 
P 



,P = P T >0,Pe: 



Then, 



Afsin(0) —Aj cos (0) 
Af sin(0) 



A- cos(0) 

P 
P 





p 0" 




P 



(AjP 



K-AjP 



Asin(0) Acos(0) 
— Acos(0) Asin(0) 

-PAf)sin(0) (AjP — PA,) cos(0) 
f PA,-) cos(0) (AfP + PA,) sin(0) 

Vi = l,2 



<0, 



<0, 



Therefore, the theorem is proved. 

Theorem 10. Consider two stable fractional order systems 
D a x — A\x and D a x = A 2 x with order a, < a < 1, then the 
following statements are equivalent: 

(1) |arg(det(M-jW))-arg(det(^ 2 -jW))| < f V 00. 

(2) A 1 and A 2 are quadratically stable, which means that BP = 
P T >0eM" xn such that 



-(-A,)z^) r p- 



P _(_a,.)3=s <0,V/=1,2 



where = — (—A,) 2 - a , V; = 1,2 and / is the identity matrix. 

Proof 4. Define c,-(s) = det(ja^ — si), i = 1,2. According to 
Theorem 6 and common quadratic stability theorem for frac- 
tional order system with order a, < a < 1, i.e., Theorem 8, 
proof is straightforward. 

Although the theory developed in the frequency domain no 
necessarily proves the SPRness, a relation was obtained as an 
equivalent issue of quadratic stability. Concerning the ease of 
designing fractional order controllers in frequency domain, the 
stability analysis in frequency domain will be really useful for 
fractional order switching systems. 

4. ILLUSTRATIVE EXAMPLES 

In this section, some examples are given in order to show 
the applicability of the theories developed for fractional order 
switching systems. 

Example 1. Let us consider the switching system (8) with 

0.3529 1.6044 



order a = 0.6, where Ai = 



0.3661 0.9237 
-0.4618 -0.1558 



1.6044 -4.4602 
Applying Theorem 8 yields: 



and A 2 = 



1 

-1 -3 





-0.5 



1 

-0.5 



Then, choosing a common matrix P 
conditions 

T 



3 1 
1 4 



, the stability 



(-MY 



p+p 



•(-Ai) 



I 

2-a 



-2 
-4 



-4 
-22 



.(-A 2 )**) P + p 



(-A 2 ) 2 -« 



-1 0.5 
0.5 -2 



<0 



<0 



are satisfied and the switching system is quadratically stable. 
Now let us compare the results with the frequency domain 
analysis. Applying Theorem 10, the following condition 



tan 



-1 



3co 



1 



tan 



-1 



0.5co 
0.5 + co 2 



<?,Vo) 
2 



(27) 



should be satisfied. The phase difference of (27) is depicted 
in Fig. 1. As can be observed, the maximum phase difference 
is 51.51°, which is less than 90°, that implies the switching 
stability condition is satisfied and the system is quadratically 
stable. 




10" 



10" 10" 

Frequency (rad/s) 
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Fig. 1. Phase difference of condition (27) for the system in 
Example 1 



Example 2. Now, let us consider the switching system given 
by (8) with order a = 0.75, where A\ = " ' 



1 
-0.5 0.1 

sf\ = 



-1 0.5 

Applying Theorem 10, we have: 



and An 



-0.3684 1.0263 
-1.0263 0.1448 



-0.2450 1.0390 
-0.5195 -0.1411 



and the following frequency domain condition 
'0.2236co\ , / 0.3861a) 



tan 



l+O) 2 J 



tan 



0.4977 + co 2 



< -,Vo) (28) 
2 



should be satisfied. The same as previous example, condition 
(28) is depicted in Fig. 2. It is shown that the maximum phase 
difference of condition (28) is 80.02°, so the switching system 
is quadratically stable. 

Example 3. Let us now consider the same switching system as 
in Example 2, but with an order bigger than 1, 1 < a < 2. 
Applying Theorem 9, the following condition: 




(29) 



where (j) = should be satisfied for all a, 1 < a < 2. Figure 3 
represents the condition (29) when the fractional order a is 
changing in the interval (1,2). In order to make this example 
clearer, the interval of variation of a is divided into three 
subintervals. As a matter of fact, Fig. 3 (a) shows the phase 
difference (29) for systems with the order a € (1, 1.5] U [1.7,2), 
whereas Fig. 3 (b) corresponds to systems with order a 6 
(1.5,1.7). As can be seen, the system is quadratically stable 
if its order a E (1,1-5] U [1.7,2). The stability region of the 
considered system is shown in Fig. 4, in which the maximum 
values of (29) are plotted versus its order a. 

5. CONCLUSION 

This paper studies the quadratic stability for fractional order 
switching systems. In particular, equivalent Lyapunov condi- 
tions in frequency domain are developed for this kind of sys- 
tems to prove their quadratic stability. Some illustrative exam- 
ples are given to show the applicability and validation of the 
proposed theory. 

Our future efforts will focus on finding a relation between the 
frequency domain method proposed in this paper and SPRness. 
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